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Abstract 

We have used Minor Planet Center (MPC) data and tools to explore the discovery circumstances and properties of the 
currently known population of over 10,000 NBAs, and to quantify the challenges for follow-up from ground-based 
optical telescopes. The increasing rate of discovery has grown to ~ 1,000/year as surveys have become more sensitive, 
by 1 mag every ~7.5 years. However, discoveries of large {H < 22) NBAs have remained stable at ~365/year over 
the past decade, at which rate the 2005 Congressional mandate to find 90% of 140 m NBAs will not be met before 
2030 (at least a decade late). Meanwhile, characterization is falling farther behind: Bewer than 10% of NBAs are well 
characterized in terms of size, rotation periods, and spectral composition, and at the current rates of follow-up it will 
take about a century to determine them even for the known population. Over 60% of NBAs have an orbital uncertainty 
parameter, f/ > 4, making reacquisition more than a year following discovery difficult; for H > 22 this fraction is 
over 90%. We argue that rapid follow-up will be essential to characterize newly-discovered NBAs. Most new NBAs 
are found within 0.5 mag of their peak brightness and fade quickly, typically by 0.5/3.5/5 mag after 1/4/6 weeks. 
About 80% have synodic periods of <3 years that would bring them close to Barth several times per decade. However 
follow-up observations on subsequent apparitions will be difficult or impossible for the bulk of new discoveries, as 
these will be smaller (H > 22) NBAs that tend to return lOOx fainter. We show that for characterization to keep pace 
with discovery would require: Quick (within days) visible spectroscopy with a dedicated >2 m telescope; long-arc 
(months) astrometry, to be used also for phase curves, with a >4 m telescope; and fast-cadence (<min) light curves 
obtained rapidly (within days) with a >4 m telescope. Bor the already-known large (H < 22) NBAs, that tend to return 
to similar brightness, subsequent-apparition spectroscopy, astrometry, and photometry could be done with 1-2 m 
telescopes. 
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1. Introduction 

Near-Barth Asteroids (NBAs) are asteroids that have 
been brought into the inner Solar System mainly 
through gravitational interactions with Jupiter and Sat¬ 
urn, placing them in orbits with perihelia q < 1.3 AU 
that intersect or come close to Barth’s orbit (Shoemaker 
et al., 1979; Bottke et al., 2002; Greenstreet et al., 2012). 
This population is of interest to science as it provides a 
probe into the dynamical and compositional evolution 
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of our Solar System and, not least, as a source of ob¬ 
jects that have shaped Barth’s history, geologically and 
biologically, through numerous impacts over many mil¬ 
lions of years. And they may yet do so again (Shoe¬ 
maker, 1983). 

NBAs have attracted public interest in recent times 
due to the unexpected 2013 Bebruary 15 event when a 
~ 17 m meteor exploded with 500 kt TNT of energy over 
the city of Chelyabinsk, Russia and injured over 1,000 
people (Bmel’yanenko et al., 2013; Borovicka et al., 
2013). The start-up of two private companies aiming 
to mine valuable resources from NBAs in the near fu¬ 
ture has rekindled interest in the scientific and for-profit 
exploration of asteroids (e.g., Blvis, 2014). In addition, 
the 44^^ U.S. president, Barack Obama, has made NBAs 
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the prime targets for human space exploration^. 

In 2005 the U.S. Congress issued a mandate to NASA 
to find at least 90% of the NBAs larger than 140 m 
(corresponding to an absolute magnitude^ ofH< 22) 
by 2020 (Brown Jr., 2005). There are an estimated 
13,200 + 1,900 in this size range, and 20,500 + 2,000 
larger than 100 m, or ff < 23 (Mainzer et al., 2011). 
Some argue that follow-up efforts aimed at characteriza¬ 
tion should concentrate only on the subset of NBAs that 
come close to Barth (Minimum Orbit Intersection Dis¬ 
tance^, MOID < 0.05 AU) and are larger than -140 m; 
these are formally classified as Potentially Hazardous 
Asteroids (PHAs)^. However, there are good reasons to 
design follow-up programs that are broader than that. 
In terms of hazard assessment, we know that the much 
more numerous objects smaller than 140 m can pose 
significant risks (Brown et al., 2013), as we were re¬ 
minded by the Chelyabinsk event. A large majority 
(-67%) of the >4,900 known NBAs with H > 22 have 
MOID < 0.05 AU. In addition, because PHAs (irre¬ 
spective of size) are, as far as we know, drawn from the 
same general population as other NBAs, characterizing 
the general population will allow us to calibrate the rela¬ 
tionships needed to infer physical and orbital evolution 
properties of the PHA subset. Binally, from the techno¬ 
logical standpoint, NASA’s proposed Asteroid Redirect 
Mission^ requires an 8 m NBA in a favorable orbit be 
found. A second option for ARM is the retrieval of a 
1-5 m boulder from a larger NBA (Abell et al., 2014). 

As we argue below, most kinds of follow-up that 
would help to quantify an object’s threat or usefulness 
must be done within a short window of time, before 
some of the properties are well-enough known to judge 
whether it is an “important” object for follow-up or not. 
Thus we concentrate our discussion on the NBA popu¬ 
lation as a whole. 

In this paper, we use the discovery circumstances and 
properties of the currently known population of NBAs 
to quantify some of the challenges that ground-based 
observers face in making follow-up observations. Our 


^http://WWW.nasa.gov/news/media/trans/obama_ksc_ 
trans.html 

^for asteroids, the absolute magnitude H is given by the apparent 
V magnitude that the asteroid would have if it could be observed from 
1 AU away, at zero phase angle, while it was 1 AU from the Sun. 

^The Earth MOID is the shortest distance separating the orbit of an 
asteroid from that of Earth. A small MOID value does not necessarily 
imply a risk of impact as both Earth and the asteroid are rarely at the 
points of their orbits closest to each other at the same time. 

^http://www.minorplanetcenter.net/iau/lists/ 
Dangerous.html 

^http://www.kiss.caltech.edu/study/asteroid/ 
asteroid_final_report.pdf 


goal is to define the optimal follow-up strategies for 
NBAs that would allow the bulk of the discovered pop¬ 
ulation to be characterized on a one-decade timescale. 
This is desireable from a planetary defence perspective 
because the bulk of the impact hazard resides with the 
smaller {H < 22), and thus more numerous, objects, 
which are also the least characterized. Gathering these 
data (composition, structure, size, etc.) for a significant 
proportion of this size population would allow for better 
estimates of potential damage due to impact, and more 
optimistically, better designed missions to defiect or de¬ 
stroy them. Bor example, if we knew that most 25 m 
NBAs are monolithic, we would produce different im¬ 
pact damage estimates and deflection mission require¬ 
ments than if we discovered most of them are rubble 
piles. Scientifically, knowing the spectral class and spin 
distributions of a significant portion of a given size pop¬ 
ulation will enable tests of current dynamical evolution 
models of the inner Solar System as well as possible 
insights into the collisional history of asteroids and the 
effect of Solar radiation pressure on spin axis evolution. 

We describe our data sources in §2, and present cur¬ 
rent discovery and characterization trends in §3. We 
then describe the known NBA population in terms 
brightness behavior (§4), positional uncertainties (§5), 
sky motions (§6), and hemisphere bias (§7). In §8, we 
discuss the implications for follow-up astrometric, pho¬ 
tometric, and spectroscopic observations. Binally, in §9 
we summarize by laying out strategies for increasing the 
rate of characterization. 

In this paper, we use the discovery circumstances and 
properties of the currently known population of NBAs 
to quantify some of the challenges that ground-based 
observers face in performing follow-up observations. 
Our goals are to define the optimal follow-up strategies 
for NBAs so the bulk of the discovered population can 
be characterized on a one-decade timescale, and to pro¬ 
vide a centralized source of NBA discovery properties 
useful for observers planning follow-up programs. 

2. Method 

The lAU Minor Planet Center^ (MPC) serves as the 
world’s clearinghouse for NBA observations and or¬ 
bital data, and provides web-based tools for generating 
lists of objects that are currently observable, and pro¬ 
viding NBA ephemerides. We have written a series of 
Python programs making use of modules Mechanize^ 


^http://WWW.minorplanetcenter.net 
^http://wwwsearch.sourceforge.net/mechanize/ 
download.html 
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and BeautifulSoup4^, to extract data from the MPC 
webpages and files. We used the following MPC re¬ 
sources. 

Data: 

• NEA. txt^ provides orbital elements at a nominal 
epoch, along with absolute magnitude //, slope pa¬ 
rameter number of observations, number of 
oppositions, arc length (for single opposition ob¬ 
jects) or years of first and last observations (for 
multiple opposition objects), and orbital uncer¬ 
tainty parameter U (see §5.1 for an explanation of 
this parameter). 

• The lists of Atens^\ Apollos^^, and Amors^^ pro¬ 
vide much of the same data but also include dis¬ 
covery date and site, and Earth MOID, but exclude 

U. 

Tools: 

• MPEph^^ is used to generate a table of ephemerides 
for user-selected NBAs for a range of dates and 
times. The data returned include V magnitude es¬ 
timates. 

• NEAobs^^ is used to generate a list of NBAs suit¬ 
able for observation according to user-specified 
criteria including observing location, NEA magni¬ 
tude, sky motion, solar elongation, RA, Dec, and 
orbital uncertainty. 

We have also made use of the web services of 
NEODyS-2^^, which provides convenient text tables 
of observations for each NEA taken from the MPC’s 
database. 

To quantify the various observational follow-up chal¬ 
lenges, we use these resources to look at the proper¬ 
ties of all NBAs discovered as of 2013 March as de¬ 
scribed in the sections below. Eor some of the analysis 


^http://pypi.python.org/pypi/beautifulsoup4/ 

^http:// WWW .minorplanetcenter.net/iau/MPCORB/NEA. 
txt 

^^The Slope Parameter, G, is a measure of how an object’s bright¬ 
ness surges when it nears opposition, the so-called opposition ef¬ 
fect. It is believed to be an interplay of shadowing and coherent- 
backscattering mechanisms (e.g., Muinonen et al., 2002). 

^^http ://WWW .minorplanetcenter.net/iau/lists/ 
Atens.html 

^^http ://WWW .minorplanetcenter.net/iau/lists/ 
Apollos.html 

^^http ://WWW .minorplanetcenter.net/iau/lists/ 

Amors.html 

^^http ://WWW .minorplanetcenter.net/iau/MPEph/ 
MPEph.html 

^^http://scully.cfa.harvard.edu/cgi-bin/neaobs.cgi 

^^http://newton.dm.unipi.it/neodys/ 


we also make use of a subsample consisting of the 6,763 
NBAs discovered during the 10 years spanning January 
1, 2002 to December 31, 2011. We call this the “Decade 
Sample.” Bor each NEA in the Decade Sample, we gen¬ 
erated a file containing a list of daily ephemerides for 
that time span, including calculated apparent V magni¬ 
tudes. 

3. Current NEA Characterization Progress Rates 
and Trends 

The current discovery rate of NBAs is ~ 1,000/year 
(Bigure 1) and this rate should increase shortly to 
> 1,500/year with the enhanced Catalina Sky Survey 
detector arrays being installed (E. Christensen, priv. 
comm.), the increased use of Pan-STARRS-1 for NEA 
discovery (100% of available time beginning March 
2014) and the coming online of Pan-STARRS-2 (R. 
Wainscoat, priv. comm.), and the first light of the BIN- 
EAR SST for NEA searches (T. Spatir, priv. comm.). 
We note that while the discovery rate has increased in 
absolute number almost continuously since 1998, this 
increase has been due mainly to an increase in the dis¬ 
coveries of the smaller H > 22 NBAs. Meanwhile, the 
discovery rate of H < 22 NBAs has remained fiat at 
~365/year thanks to a continuous decrease in the dis¬ 
covery ofH< 18 NBAs being offset by a similar in¬ 
crease in the discovery of 18 < // < 22 NBAs (see 
Bigure 1). At this rate (and bearing in mind discov¬ 
ery rate will drop as the number of undiscovered ob¬ 
jects decreases and if survey capabilities remain un¬ 
changed), the Congressional mandate of 2005 (Brown 
Jr., 2005) to find at least 90% of the NBAs larger than 
140 m (H < 22, an estimated 13,200 + 1,900 accord¬ 
ing to Mainzer et al., 2011) will be met no earlier than 
2030, i.e., more than 10 years late. Mainzer et al. (2011) 
also predict there are 20,500 + 2,000 NBAs larger than 
100 m{H< 23). 

While the pace of NEA discoveries has increased 
each year, characterization is falling farther behind. 

3.1. Sizes 

Only a small fraction (~7%) of NBAs have well- 
determined diameters and albedos from thermal IR de¬ 
tection with NEOWISE (Mainzer et al., 2011) and 
Spitzer (Thomas et al., 2011) and fewer than 5% have 
been detected by ground-based radar^^. While a higher 
radar rate (~ 100/year) has now been implemented and 


^^http://echo•jpl•nasa.gov/asteroids/index.html, ac¬ 
cessed 2014/05/01. 
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Year 

Figure 1: Discovery rate of NBAs for 1994-2013. Blue 
denotes H <22 and yellow, H > 22. 


higher spatial resolution is being investigated, the total 
number of NBAs accessible to radar is limited (by the 
outgoing pulse power) to those that come well within 
~0.1 AU*l Although -61% of NBAs have MOID < 
0.1 AU, the correct positioning to achieve a <0.1 AU 
distance from Barth occurs rarely. With the revived NB- 
OWISB mission, the prospect of more thermal IR de¬ 
tections has improved. In 3 years -2,000 known NBOs 
(Near-Barth Objects) will be detected and then albe¬ 
dos measured to -50%, for a total of -25% of known 
NEOs'^ 

Bor the vast majority of NBAs these quantities must 
be inferred from ground-based visible/near-IR observa¬ 
tions yielding estimates of absolute magnitude H and 
albedo. However, the H alone is often uncertain by 
up to a magnitude due to inhomogeneous photome¬ 
try and star catalog magnitude biases (Williams, 2015), 
light curve variability (amplitudes typically -0.4 mag, 
Mainzer et al., 2011), poorly-constrained phase curves, 
and systematic size-dependent biases of unknown cause 
(e.g., Chesley et al., 2002; Pravec et al., 2012). 

3.2. Rotation Periods 

By the beginning of 2014, <8% of all known NBAs 
had period determinations listed in the Bight Curve 


^^See the SNR graphs for the Arecibo (http://echo.jpl. 
nasa.gov/~lance/snr/far_asnrl8.gif) and Goldstone (http: 
//echo.jpl.nasa.gov/~lance/snr/far_gsnr-20.gif) radar 
dishes 

^^http://www.Ipi.usra.edu/sbag/meetings/j an2014/ 
presentations/08_1415_Mainzer_SBAG.pdf 


Database (Warner et al., 2009), with the number increas¬ 
ing at a rate of 150-200 NBAs per year (B. Warner, 
priv. comm.), mostly for V < 18 objects. We note that 
many of the light curves these rotation periods are de¬ 
rived from are not published (e.g., the MFC receives 
light curve data for only -30 NBAs per year^^). 

3.3. Composition 

Spectral characterization of NBAs proceeds at a sim¬ 
ilarly slow rate. The largest program of optical-near- 
IR spectroscopy (0.5-2.5 microns) is the MIT-IRTB 
programme^\ which has acquired a total of -1,000 
spectra (of mostly large if < 15 NBAs) at a rate of 
-100/year. It would take almost two centuries to com¬ 
plete spectroscopy of all -20,500 NBAs with if < 23 
and -120 years for the -13,200 NBAs with if < 22 
(Mainzer et al., 2011) at current rates, should these be 
the goals of such an IR spectroscopy program. It is 
noted, as we will show in §4.2, that most of these NBAs 
will not be bright enough for ground-based IR spec¬ 
troscopy. 

3.4. Orbital Uncertainties 

Orbital uncertainties are large (greater than 6 ar- 
cmin/decade along the orbit) for nearly all objects 
smaller than H = 22 owing to a lack of long-arc astrom¬ 
etry. And, as we show in §5.1, astrometric follow-up is 
also lagging. 

3.5. Trends 

Bollow-up observations will become even more chal¬ 
lenging as surveys become sensitive to fainter objects 
and as large NBAs become fully known, as they ef¬ 
fectively already are for objects with diameters >1 km 
(Mainzer et al., 2011). The progression of NBA discov¬ 
ery magnitudes (see §4) is shown in Bigures 2 and 3. 
Bor the past three decades there has been a steady in¬ 
crease in mean discovery depth by about 1 magnitude 
every -7.5 years, from (V) = 16 to (V) = 20. The mean 
H magnitude has become correspondingly larger, from 
(ff) - 15 in 1980 to (if) - 22 in 2013, corresponding 
to an approximate size decrease from 3,550 m to 140 m. 

Bigure 4 shows in more detail the breakdown of dis¬ 
covery magnitudes over the past 20 years. The discov¬ 
ery rate at V < 18 saturated at -200/year in 2000 as 
most of the big objects were found. Meanwhile, the dis¬ 
covery rate at V = 20-21 has increased over the same 
period and now dominates at -375/year. 


^*^http: //www. minorplanet center. net/light_curve 
^^http://smass.mit.edu/minus.html 
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Figure 2: Mean discovery magnitude for NBAs (ap¬ 
proximately V band, see §4.1). Bars represent l-cr scat¬ 
ter around the mean. There is a steady 1 mag increase in 
depth every ~7.5 years visible in the past few decades. 
The current mean is 20.0 with l-cr scatter of 1.0 mag. 



Figure 3: Mean H magnitude of discovered NBAs by 
year. Bars represent l-cr scatter around the mean. 


4. NEA Brightnesses 

4.1. Brightness at Discovery 

To characterize the brightnesses of NBAs at discov¬ 
ery, we have combined the MFC Atens/Apollos/Amors 
files, extracted the discovery dates of the nearly 10,000 
NBAs known as of 2013 March, and looked up their 
discovery observation optical magnitude through the 
NBODyS-2 interface. These are usually V or R magni¬ 
tudes, which we use with no correction given our coarse 
magnitude binning, but we note that we might expect 



Bigure 4: Discovery magnitude (approximately V band, 
see §4.1) for objects discovered since 1994. Lines are 
labeled by magnitude m representing the range m < 
Vdisc < w + 1. The changing slopes indicate how 
the brighter discoveries have been saturated while the 
fainter ones are increasing. The annual counts are use¬ 
ful for planning follow-up programs. 


y - /? as big as ~0.4-0.5 mag based on Sloan Digital 
Sky Survey (SDSS) color-color plots for main belt ob¬ 
jects in Ivezic et al. (2001) and transformation to R and 
V through the SDSS photometric relations^^. We also 
note that Pan-STARRS discoveries have been reported 
in, e.g., the r, /, and wide w bands, but the offsets from 
these to V are expected to be 0.1-0.4 mag for most as¬ 
teroids (Denneau et al., 2013). There are a few objects 
(currently <0.2%) that were WISE discoveries with no 
magnitude for the discovery date. Otherwise we expect 
that our technique provides a reasonable approximation 
of V band magnitude at discovery for nearly all the ob¬ 
jects. The discovery magnitude distribution is shown in 
Bigure 5 for various H bins. The overall mean is 19.4 
with a standard deviation of 1.4. 

To explore how a typical NBAs brightness changes 
during its discovery apparition, we consider the Decade 
Sample. Bor each NEA in this sample, we look up the 
discovery date in the Atens, Apollos, and Amors data 
files, and then examine the ephemerides at daily inter¬ 
vals to determine how the asteroid’s magnitude (theoret¬ 
ically) changed throughout its discovery apparition. We 
define Vdiff = ilVdisc-Iminl, where Vdisc is the discovery 
magnitude and Vniin is the brightest magnitude achieved 
during the apparition, and where -r/- values represent 


^^http://www.sdss.org/dr7/algorithms/ 
sdssUBVRITransf orm.html 
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Figure 5: Discovery Magnitude. First recorded optical 
magnitude (which we use as a proxy for V ; see text) of 
NBAs known as of March 2013 (9,723 objects). The 
current mean discovery magnitude, V ~ 20, is shown 
for reference. 


NBAs brighter after/before discovery, respectively. Fig¬ 
ure 6 gives the distribution of Vdiff. NBAs discovered 
before reaching peak brightness have positive values of 
VdifF, while those discovered after have negative values. 
Most NBAs (~60%) are discovered within 0.5 magni¬ 
tudes of their peak brightness during the discovery ap¬ 
parition. 



Figure 6: Distribution of the difference between discov¬ 
ery and brightest magnitude, Vdiff, for the 6,763 NBAs 
in the Decade Sample. NBAs discovered before reach¬ 
ing peak brightness have positive values of while 
those discovered after have negative values. 57 NBAs 
with I Vdiff I > 4 were omitted for clarity. 


We then evaluated the fraction of NBAs remaining 
brighter than a series of limiting apparent magnitudes, 
Viim, after discovery (Figure 7). It is clear that this frac¬ 
tion drops rapidly for many interesting cases. The val¬ 
ues of Viim = 18 was motivated by the infrared spec¬ 
troscopic capabilities of the IRTF, and also is in range 
for photometric and astrometric follow-up by amateurs. 
Viim = 20 represents the infrared spectroscopic limit 
of Magellan (see §8) and is the current mean discov¬ 
ery magnitude. Vum = 21 represents the limit for the 
proposed NBA spectroscopic follow-up project LIN- 
NEAUS (Elvis et al., 2014). The solid lines show 
the brightness evolution using the full Decade Sample, 
while the dashed lines are drawn from only the NBAs 
discovered during the last 2 years of the Decade Sam¬ 
ple. As can be seen, the number of NBAs that remained 
above a certain brightness decreased noticeably over 
time, due mainly to the increase in the number of small 
(H > 22) NBAs discovered (Bigure 1) and the greater 
overall survey depth in more recent years (Bigure 2). 



Days 


Bigure 7: Observable days following discovery for the 
6,763 NBAs in the Decade Sample. Percentage of 
NBAs that remained brighter than V = 18, 19, 20, 21, 
22, 23, 24 and 25 for up to 10 weeks following their 
discovery. Solid lines are for the full Decade Sample; 
dashed lines are for the last 2 years of the sample only. 


Approximately 2/3 of the objects discovered during 
the decade were found at V < 20, but faded by 0.5 mag 
after ~1 week, 3.5 mag by ~4 weeks, and by more 
than 5 mag by ~6 weeks. Meanwhile, at the bright 
end, we see that only -17% of NBAs became brighter 
than V = 18, even during their discovery apparition 
(Bigure 7), and even less, -3%, remained brighter than 
V = 18 for more than 30 days. The small objects 
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{H » 22) that are only seen on very close passages 
would fade much faster than these aggregate values, as 
evidenced by the dashed lines, which show the values of 
these curves for the last 2 years of the Decade Sample. 
These last 2 years featured many more small {H > 22) 
NBA discoveries than for the average of the sample. 
The discovery apparition then offers only a short win¬ 
dow for characterization. We explore the implications 
for follow-up in §8. 


4.2. Brightnesses on Subsequent Apparitions 

NBAs will have other apparitions following discov¬ 
ery, so we also investigated their observability at these 
subsequent apparitions. The success of a “subsequent 
apparition” approach to characterization depends on 
two factors: 


1. The length of time that will elapse before the 
NBAs next apparition. 

2. The brightness the NBA will reach on its next ap- 
parition(s). 

The synodic period, 5, is the time that elapses be¬ 
tween two consecutive apparitions (i.e., when the NBA 
returns to a similar position near Barth). If an NBAs 
orbital period is denoted by P (given in years), then its 
synodic period with respect to Barth will be given by 
(Breedman and Kaufmann, 2005): 


S = 



years 


( 1 ) 


The distribution of synodic periods for the known 
NBAs as of March 2013 is shown in Bigure 8. 

Almost 80% of known NBAs have S < 3 years 
and will therefore have at least three apparitions each 
decade. Synodic periods of 3-10 years account for 
~20% of known NBAs. These reapparition times are 
too long for space mission target selection studies as the 
properties and highly accurate orbits needed for selec¬ 
tion will take too long for planning purposes. 

As an example of what it takes to determine an or¬ 
bit well enough for space mission planning, we present 
NBA (101955) Bennu, the target for the OSIRIS-RBx 
mission^^. To obtain an accurate enough ephemeris, it 
was observed with radar during 3 apparitions and had 
optical observations spanning 15 years. This provided 
a very accurate orbit determination, which also con¬ 
strained the magnitude of the Yarkovsky effect acting 
upon it, thus allowing the mass/area ratio of Bennu to 
be calculated (Chesley et al., 2014). 


^^http://WWW.asteroidmission.org/ 



Synodic period (years) 


Bigure 8: Synodic Period distribution of NBAs known 
as of March 2013 (9,723 objects). Bor clarity, the 5% of 
outliers with S > 10 years are not plotted. 


Only 5% of the nearly 10,000 NBAs known as of 
2013 March have S > 10 years as their position relative 
to Barth changes very slowly; it will be decades before 
they can be re-observed from Barth after their discov¬ 
ery apparition. There may also well be NBAs with long 
synodic periods currently “hiding” within the region of 
space at elongations <60°, making them unobservable 
because of their position with respect to the Sun; these 
may not be discovered from Barth for decades. The pos¬ 
itive side of this situation is that once a long synodic pe¬ 
riod NBA becomes observable, it will remain so for an 
extended period of time. 

Given that most of objects will return several times 
within a decade, we computed the brightest magnitude 
that NBAs known as of 2013 March will achieve on re¬ 
turn apparitions over a 10 year period starting January 
1, 2018, by which time, most of today’s newly discov¬ 
ered NBAs would be making at least their first return. 
We used the MPC tool MPEph (see §2) to generate daily 
ephemerides for these NBAs and determined the bright¬ 
est V magnitude reached over a 1, 5, and 10 year period 
(starting January 1, 2018). Results are shown in Big¬ 
ure 9 (on the left for objects with H < 22 and on the 
right for H > 22). 

Bigure 9 shows that over a period of one year, less 
than 1% of the known NBA population will become 
brighter than y = 18. Bven over a decade, only -17% 
will reach V < 18. By comparison to Bigure 5, we 
see the larger NBAs (H < 22) typically do return to 
similar brightnesses within a decade; however, smaller 
NBAs (H > 22) are generally no brighter than 5 magni- 


7 




Figure 9: Follow-up apparitions of known (as of March 2013) NBAs over a 1, 5 and 10 year period (starting January 
1, 2018) for (left) H < 22, and (right) H > 22, showing the distribution of the brightest V magnitude that will be 
reached. The current mean discovery magnitude, V~20, is shown for reference. Note that fainter objects are to the 
right in each plot. 


tudes fainter than at discovery, and rarely brighter than 
y = 21. This reflects the bias towards discovering small 
NBAs during their very closest approaches to Barth. 

Most of the NBAs that remain to be discovered and 
characterized are at the smaller end of the size distri¬ 
bution^^; if follow-up of these objects is not performed 
during the discovery apparition, the distribution in Big- 
ure 9 suggests that nearly all will remain uncharacter¬ 
ized for over a decade. We conclude that from a bright¬ 
ness point-of-view, follow-up cannot wait for a subse¬ 
quent apparition but must be performed during the dis¬ 
covery apparition. 

5. NEA Positional Constraints 

5.7. Orbit Uncertainty 

Bollowing discovery, an NBAs orbit must be deter¬ 
mined to sufficient accuracy for the NBA to be success¬ 
fully re-acquired for follow-up on a subsequent appari¬ 
tion. The uncertainty in an NBAs orbital solution is 
quantified by the MFC via the Orbit Uncertainty Param¬ 
eter, . U is an integer running from 0 (extremely 
precise orbit) to 9 (highly imprecise orbit). An asteroid 


^^http://sservi.nasa.gov/wp-content/uploads/2014/ 
03/Harris.pdf 

^^http://www.minorplanetcenter.net/iau/info/ 
UValue.html 


with a U value of 0 is defined to have an orbit precise 
enough to allow calculation of its position along its or¬ 
bital path to within +1" after one decade. Bach integer 
step in U is a logarithmic step increase of a factor ~4.4 
in this uncertainty. So U = 9 corresponds to an uncer¬ 
tainty of more than -41° after a decade. It is important 
to understand that U is associated with an uncertainty in 
position along the orbit of the asteroid, which cannot be 
readily converted to an on-sky positional uncertainty be¬ 
cause it would depend on the asteroid’s trajectory with 
respect to the Barth observer at the time of the observa¬ 
tion. 

With sufficient observing time and telescope aperture, 
any NBA of sufficient brightness can be recovered at 
subsequent apparitions after discovery, but practically 
speaking, to lie within a typically-sized imaging field 
(-10 arcmin), requires U < 3 (for an NBA; for Main 
Belt Asteroids the tolerances are more relaxed). This 
quality of U has to be reached during the discovery 
apparition to allow for ready recovery, and thus fur¬ 
ther orbit quality improvement at subsequent appari¬ 
tions. Radar has the potential to greatly improve the 
quality of a calculated orbit given that a radar observa¬ 
tion can reduce on-sky positional uncertainty of an NBA 
at the following apparition by up to 2 orders of magni¬ 
tude or more (Ostro et al., 2002, 2007). Unfortunately, 
only -100 NBAs are observed each year with radar and, 
as we discuss in §3.1, many NBAs are out of reach of 
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current radars. 
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Figure 10: Orbit Uncertainty vs. if magnitude for NBAs 
discovered before 2013 March. A line marking H = 22 
is plotted for convenience. 


We used the MFC NEA. txt to extract the values of U, 
if, number of oppositions, and arc length (given as days 
for single opposition objects, or years between first and 
last observation for multi-opposition objects). Figure 10 
shows the U values as a function of the if magnitude for 
NBAs discovered through 2013 March^^. Over 65% of 
larger (H < 22) NBAs have values of U <3, and thus, 
orbits well-defined enough for re-acquisition on subse¬ 
quent apparitions. These high-quality orbit determina¬ 
tions come about for two reasons: (1) many were dis¬ 
covered over a decade ago so they have been observed 
over several apparitions; and (2) being large, they are 
generally brighter and thus can be discovered (or re¬ 
observed) farther from Earth, where their relative mo¬ 
tion to Earth leads to slower fading than is the case for 
smaller objects traveling closer. 

However, for smaller NBAs (if > 22) less than 10% 
have U < 3. By fi > 25 this fraction is reduced to ~5%. 
Only more recent surveys have been able to find signif¬ 
icant numbers of H > 22 NBAs (Figure 1). In addition, 
this population tends to return at fainter magnitudes (see 
§4.2 and Figure 9). Therefore, few have been observed 


^^The estimated diameter was calculated using the following rela¬ 
tionship (Fowler and Chillemi, 1992) between the diameter, D, visual 
albedo, p, and absolute magnitude, H: 

1329 

D = -10“®'^^ kilometers (2) 

Vp 

We used a value of p = 0.14, as recommended by Stuart and Binzel 
(2004). 


over multiple apparitions. The overwhelming fraction 
of known H > 22 are thus poor candidates for targeted 
re-acquisition on the next apparition and will have to 
be re-discovered by chance. There are only ~50 very 
small NBAs {H > 25) with very well-determined orbits 
(U < 2) such that they could be found beyond the next 
apparition with a targeted observation. 

At if > 26 the fraction of NBAs with smaller U 
shows a small increase. This we attribute to the fact 
that smaller NBAs must be physically close to Earth to 
become bright enough to be detected. When an NEA is 
discovered close to Earth, its initial orbit will not gen¬ 
erally be accurate enough to entirely rule out an im¬ 
pact with our planet. Possible impactors are heavily 
observed until an improved and more accurate orbit is 
obtained, with a sufficiently low value of U to rule out 
an impact with 100% certainty. 

5.2. Arc Length 

The most direct way to reduce the uncertainty in 
an NEA orbit is to increase the arc length observed 
during the discovery apparition, and then through re¬ 
observation at subsequent apparitions. An NBAs arc 
length is simply the number of days elapsed between 
the first and last time it was observed. 

A clear relation of U with log(arc length) in days 
is shown in Figure 11. This figure can be compared 
to Figure 4 in Desmars et al. (2013), who carried out 
a similar exercise using the current ephemeris uncer¬ 
tainty as a measure of the orbital uncertainty for all 
known asteroids (not just NBAs). They found that as¬ 
trometry conducted within the first 2 days after discov¬ 
ery had the highest impact in reducing positional un¬ 
certainty for NBAs; astrometry from that point on un¬ 
til 7-14 days had little effect on the uncertainty, after 
which point the uncertainty begins to decrease with in¬ 
creasing arc length. Our results corroborate this for the 
NBAs: Up to ~10 days, the U is typically ~7. Des¬ 
mars et al. (2013) also found that astrometry on 10-250 
day timescales improves the orbit much more rapidly. 
Our plot bears these results out for NBAs as well: U 
decreases linearly as a function of log(arc length) for 
arc lengths >10 days. There is a clear ridge line in our 
figure that defines the typical arc length needed to reach 
a given U value, and by 100-200 days the U has typ¬ 
ically dropped to ~3. The relation between U and arc 
length extends to multiple apparitions. However, even 
with multiple apparition observations, about -40% of 
NBAs still have U >2\ this is generally due to an NBA 
having been observed for short arcs during a few (of¬ 
ten non consecutive) apparitions, stressing the need for 
longer individual observing arcs during each apparition. 
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Achieving t/ < 3 at the discovery apparition typically 
requires astrometry over an arc length of >3 months 
(Figure 11). For H < 22 NBAs, obtaining such a long 
arc requires large telescopes (4 m class or larger) to en¬ 
sure the object is detected (see §8). However, the ex¬ 
tent to which further astrometric measurements reduce 
U for specific NBAs will depend on the phase of the 
orbit being observed, as determining orbit curvature is 
crucial. We will explore this in a subsequent paper. We 
did search for differences in the arc length vs U relation 
for Aten, Apollo and Amor NBAs separately, but found 
no difference. 
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Figure 11: f/ vs. arc length for NBAs discovered as 
of 2013 March. Data point density is represented by 
the color scale. Arc lengths for multiple apparition ob¬ 
jects are plotted as 365 days times the number of years 
between first and last observation. The gap at 1 year 
marks the boundary between single opposition and mul¬ 
tiple apparition objects. 


6. Sky Motion 

The relative proximity of NBAs means they have 
proper motions that can potentially limit follow-up ob¬ 
servations. In this section we explore NBA sky motions 
during a typical year. Using NEAobs (see §2) we cre¬ 
ated ephemerides for all NBAs visible during 2013, and 
determined the sky motion on the night when each NBA 
was at its brightest. Figure 12 shows the resulting dis¬ 
tribution. 

We find that the typical sky motions are <5"/min, 
with median ~1.5'7min. This plot can be used to es¬ 
timate the fraction of targets accessible to a telescope 


Figure 12: Sky Motion of NBAs visible during 2013 
(with elongations >60°). For each NBA, the sky motion 
on the night of maximum brightness was used. Outliers 
with motions >30'7min (~1% of the sample) are not 
shown. 


conducting follow-up observations from a sky motion 
perspective. For example, less than 5% of NBAs vis¬ 
ible in 2013 were moving faster than the 9"/min that 
the Magellan telescope is capable of tracking, so sky 
motions would not pose a significant limitation for this 
case. On the other hand, for any follow-up observations 
done by sidereal tracking, the sky motion constrains ex¬ 
posure times as described in §8. 


7. Hemisphere Bias 

All large NBA ground-based surveys, except the re¬ 
cently cancelled Siding Spring Survey in Australia, are 
in the Northern Hemisphere, so ~30% of the most 
southern part of the celestial sphere is not searched as 
thoroughly as the rest. Our ephemerides show that 92% 
of the known NBAs will appear in both hemispheres 
within a five year period, and 98% within a ten year 
period. For H < 22, this may permit follow-up from 
either hemisphere on subsequent apparitions. However, 
given that for H > 22 the objects will return ~5 mag 
fainter (see §4.2, Figure 9), follow-up may not be possi¬ 
ble from either hemisphere. This is investigated further 
in a future paper. 

The limitations discussed in §4.1 (see Figure 9) high¬ 
light the need to complete most of the follow-up at dis¬ 
covery. For this reason, some follow-up facilities need 
to be placed at similar locations to the telescopes being 
used for discovery. As of 2012, -67% of NBAs are be¬ 
ing discovered from the U.S. Southwest region (in an 
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area around ~32°N ~110°W) by the Catalina Sky Sur¬ 
vey and associated Mt. Lemmon Survey, LINEAR^^, 
and Spacewatch. Rapid (discovery night) follow-up 
would argue for very close proximity from a common- 
weather point of view, while longer-term follow-up 
(e.g., long-duration coverage of light curves or follow¬ 
up in subsequent nights after discovery) could benefit 
from observatories distributed along longitude and at 
similar latitudes. Follow-up of discoveries from Hawaii 
(where Pan-STARRS is currently the second most pro¬ 
lific NBA discovery site) must be made in Hawaii, as 
the next high-quality observing site to the West is in the 
Canary Islands, some 15 h away. If rapid follow-up is 
desired, then this also argues for astrometry of discov¬ 
ery observations to be submitted to the MPC in near re¬ 
altime. 

8. Discussion 

Our analysis of the currently-known population of 
NBAs shows that: 

1) The discovery rate of // < 22 NBAs has remained 
steady at ~365/year over the past decade (see Fig¬ 
ure 1), suggesting it will take over 20 years to dis¬ 
cover the remaining population of H < 22 NBAs 
unless the rate of 18 < // < 22 NEA discovery in¬ 
creases dramatically. 

2) The mean magnitude for NBAs being discovered to¬ 
day is H ~ 23, which has been increasing at a rate 
of ~1 mag every 3 years. Meanwhile, the mean dis¬ 
covery magnitude for NBAs is V ~ 20, which has 
increased at a rate of 1 mag every ~7.5 years. 

3) -60% are discovered within 0.5 mag of their peak 
brightness during the discovery apparition and they 
fade quickly: the typical NEO will dim by 0.5 mag 
after -1 week, 3.5 mag by -4 weeks, and by more 
than 5 mag by -6 weeks. 

4) 80% of synodic periods (which define windows for 
followup in subsequent apparitions) are <3 years so 
that 3 apparitions per decade are available for char- 
actrization if the returning NEA is sufficiently bright. 

5) The average brightness for H > 22 NBAs on return 
apparitions in the decade following their discovery 
is 5 magnitudes (i.e., 100 times) fainter than at dis¬ 
covery. For larger {H < 22) NBAs, about half of 
the decade-scale returns achieve similar brightness 
to that of their discovery apparition. 


^^Lincoln Lab’s Near Earth Asteroid Research program, which was 
terminated as of 2013. 


6) Orbital uncertainties are significant for the bulk of 
currently known NBAs, as most have been ob¬ 
served for relatively short arcs on a single appari¬ 
tion. Whereas most (60%) of H < 22 NBAs have 
f/ < 3, less than 10% of H > 22 NBAs do. To 
reduce the orbital uncertainty parameter on the dis¬ 
covery apparition to levels practical for recovery on 
subsequent apparitions (f/ < 3) requires astrometric 
measurements over arcs > 3 months. 

7) Sky motions when NBAs are at their observ¬ 
able brightest are generally <5'7iiiin, with median 
~1.5'7niin except for very close approaches for 
which the rate can exceed this by an order of magni¬ 
tude. 

8) Currently, the concentration of surveys in the North¬ 
ern Hemisphere means that follow-up observatories 
must be located in the North too for observations 
on the discovery apparition. However, follow-up of 
H < 22 objects on subsequent apparitions can be 
done from North or South (assuming they can be 
found), as >92% of these bigger NBAs will appear 
in both hemispheres within a 5 year period at mag¬ 
nitudes comparable to those at discovery. 

We can use these results, along with projected discov¬ 
ery circumstances for the next few years, to determine 
telescope requirements for various types of follow-up 
observations. 

8.1. Fiducial Magnitudes 

We begin by defining a few fiducial magnitudes for 
exposure time predictions. The mean magnitude at dis¬ 
covery is currently V = 20.0 +1.0 mag (see Figure 2), 
so to reach the bulk (-2/3) of these objects at discov¬ 
ery, observers would need to achieve depths of V - 21. 
However, the mean discovery magnitude among all sur¬ 
veys increases by -1 mag every 7 years, and Pan- 
STARRS is already achieving this today. Thus, we take 
y - 21 as the projected mean discovery magnitude for 
a year or two from now (and we should expect that to 
reach -2/3 of the population we would have to get a 
magnitude or so fainter, to objects with discovery mag¬ 
nitudes of y - 22). 

In the discussion that follows, we take the 
historically-preferred V magnitudes and translate them 
to R (or AB mag r), which are more likely representa¬ 
tive for the follow-up observations. We do so crudely 
by taking y - R ^ 0.5 and note that r - R ^ 0.2 (AB 
mag to Vega mag, see section §4 above). With these 
assumptions and our results from above, we derive the 
following fiducials: 
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• The fast fading of newly-discovered asteroids and 
lOOx fainter reapparitions will require some types 
of follow-up to be done rapidly (within days) af¬ 
ter discovery. Taking the projected mean discovery 
magnitude of V = 21 + 1, (/? = 20.5 + 1) and allow¬ 
ing for 0.5 mag of fading in a week, we find that 
the observations will need to achieve mean depths 
of y ~ 21.5 (R ~ 21), or for bulk of new objects 
R ~ 22. 

• For follow-up observations done on a '^months 
timescale after discovery the NBAs are expected 
to be >3 mag fainter still; i.e., we require mean 
depths of R ~ 24, or for the bulk of new objects 
R ~ 25. 

• For currently-known NEAs making return appari¬ 
tions in the next decade, we take the mean dis¬ 
covery magnitude for all NEAs discovered through 
2013 March of V ~ 19.4, with V ~ 21 encompass¬ 
ing >2/3 of the objects (see §4.1). For H < 22 
objects, the reapparitions will be similarly bright, 
setting the observing depths at V ~ 21, or 7? ~ 20.5 
for the bulk of old objects. However, for the H > 
22 objects the depths required to reach the mean is 
5 magnitudes fainter. That would ho R = 25.5 for 
the bulk of old objects, or 7? ~ 24.5 for the bright 
half of old objects. 

We show in Table 1 the signal-to-noise (S/N) and ex¬ 
posure times for magnitude ranges around our fiducial 
magnitudes under dark-sky conditions typical of moon¬ 
less nights on Kitt Peak and Cerro Tololo^^. We have 
assumed the performance of a typical modem back- 
illuminated CCD and take the radius for photometry to 
be 2 arcsec sampled with sub-arcsecond pixels. (Thus 
for excellent seeing our estimates would be pessimistic.) 
Our calculations check out against the online exposure 
time calculators^^; e.g., the Kitt Peak 2.1 m and 4 m 
telescopes with these assumptions. 

8.2. Requirements for Astrometry 

For ground-based astrometry we compute the expo¬ 
sure times needed to achieve tolerances of 0.1 arcsec. 
For faint objects, with a well-sampled PSF and no im¬ 
age trailing, astrometric precision should be given by 
Bimey et al. (2006): 


^^https://www.noao.edu/kpno/manuals/dim/, accessed 

2014/03/17. 

^^http://www.noao.edu/noao/noaonews/mar99/nodel6. 
html 


1 FWHM 

“ 2.355 ^ (S/N) ^ ’ 

In 1 arcsec seeing we then require S/N ~ 4 to achieve 
0.1 arcsec precision, but as a source at that level is 
barely significantly detected, in practice one will likely 
require a slightly higher S/N. 

The requirement for no NEA trailing could be satis¬ 
fied in several ways. One could track non-sidereally on 
the (potentially imprecisely known) motion of the aster¬ 
oid, though centroiding the resulting star trails could in¬ 
troduce subtle astrometric biases, for example through 
timing errors caused by seeing/transparency variations 
during the exposure. Alternatively, one could take sets 
of individual exposures short enough to prevent trailing 
and then combine them later, either by shift-and-add to 
optimize S/N (which the MPC sternly warns against), 
or by averaging the individual position measurements. 
For the typical NEA motion of 1.5 arcsec/min found 
above, sidereally-tracked exposures would be limited to 
r ~ 40 s. 

From Table 1 we see that a 2 m telescope could ac¬ 
complish the astrometric measurements even with side¬ 
real tracking for newly-discovered objects on timescales 
of about 1 week (not surprising since the discoveries 
will be made by 2 m-class telescopes), and for reappari¬ 
tions of known large NEAs (which will likely be done 
de facto by future search surveys). Fong-arc (months) 
astrometry of only the brightest new objects could pos¬ 
sibly be done on a 2 m by tracking on the asteroid and 
trailing the stars or shift-and-add. The mean new ob¬ 
jects would require a 4 m telescope. Fong-arc astrom¬ 
etry of the faintest new objects and the H > 22 reap¬ 
paritions would require a larger telescope still. This 
argues strongly for conducting long-arc astrometry on 
4-v m telescopes during the discovery apparitions. The 
field of view would ideally be ~20 arcmin to ensure that 
enough astrometric standards lie within each image and 
to accommodate uncertain positions. To ensure a well- 
sampled PSF would require pixel sizes of <0.5 arcsec. 

8.3. Requirements for Light Curves 

For light curves, the S/N requirements are much 
stiffer. In addition, the temporal sampling defines the 
range of rotation rates for which periods can be deter¬ 
mined, with one benchmark being that for “sparsely- 
sampled” observations over many rotations that will un¬ 
dergo Fourier analysis, the light curve needs to be sam¬ 
pled at intervals no more than 18% of the rotation period 
(Pravec et al., 2000). The nearly 700 NEAs in the Fight 
Curve Database (see Pravec and Harris, 2000; Warner 
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Table 1: S/N and exposure times for 2 m and 4 m telescopes 


R mag 

S/N 

= 4 

S/N = 

10 

S/N = 

100 

Circumstances 

2 m 

4 m 

2 m 

4 m 

2 m 

4 m 

20.5 

16 s 

2s 

50 s 

10 s 

45 min 

10 min 

Mean for new discoveries on same nights; or, 

most // < 22 reapparitions 

22.0 

90s 

16 s 

7 min 

90s 


2.5 h 

Most new discoveries at ~ 1 week 

24.0 

35 m 

10 m 

4h 

Ih 



Mean for new discoveries at ~ 1 month 

25.0 

4h 

1 h 


5h 



Most new discoveries at ~1 month; or, most 

H > 22 reapparitions 


et al., 2009) with light curve properties listed show am¬ 
plitudes of -0.4 mag, with rotation periods typically 
<2.2 h for if > 22 objects (1/4 of the sample) and >2 h 
for larger ones. 

However, the true distributions are not known, as 
NBAs with hard-to-measure periods would not be in¬ 
cluded in the database. In their attempt to measure 
light curves of 83 sub-km NBAs with a 2.4 m telescope 
(down to i? - 19), Statler et al. (2013) found that 2/3 
had ambiguous light curves that would not permit de¬ 
termination of the period. Their cadence was -90 s, 
and their median photometric precision was 0.04 mag. 
While their ambiguous cases were measured with typi¬ 
cally worse than average precision, they concluded that 
the ambiguous objects were likely either super-fast or 
super-slow rotators, or objects with intrinsically low 
amplitude (e.g., axis-symmetric) objects. 

A recent detailed study by Warner and Harris (2011) 
finds that the predicted yield of Bourier analysis on 
sparsely-sampled (up to 16 times per night) light curves 
with noise of -0.03 mag is highest in the period range 
of 2-6 h and for amplitudes >0.2 mag. Such a strat¬ 
egy could be well-suited to reapparitions of the larger 
if < 22 NBAs, typically exhibiting >2 h spin periods 
and amplitudes of -0.4 mag. Bor a 2 h period the 18% 
requirement yields spacings of >20 min. Table 1 shows 
that observations with a 1-sigma precision per point of 
a few percent, i.e., -0.03 mag (yielding light curves 
with S/N -10 and corresponding to the cases discussed 
by Warner and Harris, 2011), could be obtained with 
long blocks of time on 2 m-class telescopes. Wide-field 
surveys for transients are already providing periods for 
(mostly Main Belt) asteroids this way (e.g., Chang et al., 
2014). 

However, as Warner and Harris (2011) articulate, 
there are compelling reasons to measure the spin prop¬ 
erties of the smaller {D < 200 m) NBAs: This size 
regime has been poorly sampled; it is susceptible to ra¬ 
diation pressure effects, which could then be quantified; 


allows probing below the “spin barrier” at -2.2 h; can 
show if at these sizes, asteroids still have satellites; does 
the fraction of tumblers increase with size? 

As most new discoveries will be smaller NBAs, they 
will likely be fast-spinning and fade fast, so cannot be 
observed for many nights. Bor light curve requirements 
for these objects we adopt as a benchmark the same pho¬ 
tometric precision of 3% and set the cadence at 30 s to 
potentially reach periods under 10 minutes (or less if 
several periods are observed). We conclude that even a 
4 m-class telescope could accomplish this only for the 
brighter new discoveries on the same night. 

Some light curve work has already been done on a 
very large telescope. Kwiatkowski et al. (2010) reported 
on a pilot program to use the 10 m SABT telescope to 
measure 14 NBAs with 7? - 20, if > 21.5 with 5-60 s 
exposures (various instrumental issues limited them to 
-0.1 mag precision). They measured periods ranging 
from 77 s to 44 min. A dedicated campaign that ob¬ 
served all of the newly-discovered if > 22 NBAs for 
3 months could double the sample of periods for this 
population. 

Binally, Warner and Harris (2011) cite the potential 
usefulness (and outreach synergy) of smaller (0.5-1 m) 
telescopes in obtaining light curve observations in the 
range V = 17-19. Bxposure times to reach the re¬ 
quired S/N would force a multi-rotation Bourier analy¬ 
sis approach for most of this range on small telescopes, 
which is better-suited to the large (if < 22) NBAs. Big- 
ure 4 indicates current rates of -1 new discovery per 
night overall within this magnitude range. But for any 
newly-discovered H > 22 NBAs (and of the known pop¬ 
ulation, about 1/4 were discovered with V < 19), their 
(likely) fast rotations would demand faster cadences and 
a 2 m-class telescope. Meanwhile, for reapparitions of 
currently known NBAs, -10% (or -1,000 NBAs) will 
achieve V = 18 within a decade (§4.2). 
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8.4. Requirements for Spectra 

Optical reflectance spectra with modest resolution 
{R = 100) of Main Belt asteroids were the basis of the 
Bus and Binzel (2002) asteroid taxonomy scheme. They 
provide a breakdown of types into Carbonaceous (C), 
Stony (S), and Metallic and Unknown (X), and into sev¬ 
eral sub-types. The currently standard optical-infrared 
Bus-DeMeo taxonomy (DeMeo et al., 2009) is superior 
to the Bus-Binzel taxonomy as it adds spectral diagnos¬ 
tics in the near-IR (1.0-2.5 microns) to the optical. All 
but three of the 26 Bus-Binzel taxonomy classes are 
preserved, with just one new class (Sv) defined in the 
Bus-DeMeo scheme. For some purposes, the visible- 
only regime preserves useful discriminants among types 
(e.g., see Figures 8, 10 and 13 in Bus and Binzel, 2002). 
From the NEO perspective, many key classes are de¬ 
fined only by their visible-wavelength features; e.g., the 
selection of potential metallic M-type objects from the 
two other sub-classes of the X-class, is improved by re¬ 
jecting the -30% of E-types with a strong 0.49 mm op¬ 
tical feature. The 0.7 mm band, though uncommon in 
NEOs, is the only alternative to the 3.1 mm feature to 
diagnose hydrated minerals. 

Optical to infrared spectra of NEAs as faint as V = 

17.5 can be obtained with an exposure of 30-60 minutes 
using the spectrograph SpeX on the 3 m Infrared Tele¬ 
scope Eacility (IRTE) on Mauna Kea, Hawai’i^^ (R.R 
Binzel, priv. comm.). Eigure 4 indicates that this would 
apply to -100 (-10% and decreasing) newly-discovered 
NEAs per year. Optical to infrared spectra of NEAs 
as faint as V = 20 (Moskovitz, 2012) can be taken in 
a 1 hour exposure with the EIRE spectrograph on the 

6.5 m Magellan telescopes at Las Campanas Observa¬ 
tory in Chile^^ (Simcoe et al., 2010). However, for pro¬ 
jected new discoveries, V = 20 (R = 19.5) would reach 
only the brightest -1/3 of objects (see §8.1), and then 
only on the discovery night. Most H <22 reapparitions 
would be reachable, but only if the astrometry were of 
sufficient accuracy to be able to place the NBA in the 
spectrograph slit. 

An alternative approach is to forego infrared spec¬ 
troscopy. Optical spectroscopy is much less demand¬ 
ing than near-infrared spectroscopy as the atmospheric 
background is orders-of-magnitude fainter. A 2 m-class 
telescope with an efficient spectrograph can obtain a 
modest resolution {R = 100) spectrum with S/N -15 
for a y = 20.5 {R = 20) NBA in under 1 h (e.g., the 
LINNEAUS project; see Elvis et al., 2014); this repre¬ 
sents the mean new discovery on the discovery night, 

^^http://irtfweb.ifa.hawaii.edu/~spex/ 

^^http://obs.carnegiescience.edu/Magellan/ 


and most H <22 reapparitions. New discoveries would 
require a 4 m telescope after a week. 

This argues strongly for a dedicated 2 m telescope 
equipped with a high-efficiency visible wavelength R = 
100 spectrometer. At a good site, a telescope can ef¬ 
fectively observe faint NEAs for -210 nights of the 
year, after weather (15%), full Moon (25%) and equip¬ 
ment maintenance (5%) are excluded^^. Assuming an 
average -10 h of observation per night, a dedicated 
NBA telescope/spectrograph combination could obtain 
-2,000 spectra/year. Such a setup could keep pace 
with new NBA discoveries with same/next night spec¬ 
troscopic followup. It would double the sample of NBA 
spectra within the first 6 months and would extend the 
coverage to small objects. 

8.5. Requirements for Colors 

In the absence of spectra, broad-band colors can in 
some cases provide some level of discrimination among 
asteroid types, for example as shown for the griz SDSS 
filters by DeMeo and Carry (2013) and Ivezic et al. 
(2001). Erom their figures, we judge that photometry 
at the level of a few percent would be sufficient (higher 
precision would not necessarily be helpful as the scat¬ 
ter in colors within a population is typically on the or¬ 
der of 0.1 mag). One significant complication is that 
the NBAs light in all bands will be modulated by its 
rotation. Eor big NEAs (H < 22), periods can be 
many hours, but for the smaller NEAs periods are typ¬ 
ically <2.2 h (see §8.3). Eor very fast rotators, long 
exposures would give an average over the light curve 
and thus could be used to get (average) colors. But 
for exposures on, or separated by, timescales significant 
compared to the (generally unknown) period, there will 
be light-curve-dependent corrections depending on the 
(generally unknown) amplitude, which is often nearly 
half a magnitude (see §8.3). In any case, exposures can 
be tracked on the asteroid. 

Interpolating between the S/N =10 and 100 columns 
in Table 1, and using the R-band exposure time as a rea¬ 
sonable estimate for the average exposure time over g, 
/, and z, we conclude that 2 m telescopes could achieve 
the required precision in short (minutes per filter) ex¬ 
posures for most H < 22 reapparitions and the mean 
new discoveries. Longer exposures (hours to get multi¬ 
ple filters) would be needed for each new mean NBA 
a week after discovery, which would only be practi¬ 
cal for fast rotators but which is not likely worth the 

^^http://www.gemini.edu/sciops/statistics/ 
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time given the crude diagnostic they would permit; low- 
resolution spectroscopy as described above would be a 
better choice. A 4 m telescope would be somewhat more 
practical for multicolor photometry of fast-rotating new 
discoveries out to a week, where 4 filters could be done 
on an hour timescale. 

8.6. Requirements for Phase Curves 

By measuring the brightness of an asteroid over a 
range of phase angles (a) it is possible to derive an 
estimate of its absolute magnitude. Moreover, it has 
been shown that different taxonomic classes show dis¬ 
tinct linear relationships between phase angle and mag¬ 
nitude (for a > 3°), and that the slope of this line (the 
phase coefficient, j8) is related to the albedo of the as¬ 
teroid (Belskaya and Shevchenko, 2000). These rela¬ 
tionships appear to hold true for smaller asteroids too 
(Hergenrother et al., 2013). Thus, by measuring a phase 
curve we obtain an estimate of the asteroid’s albedo and 
broad taxonomic class, as well as its absolute magni¬ 
tude. If the absolute magnitude is already known, the 
phase curve can be used to estimate the slope parameter 
G. 

10-20% photometry would be good enough for 
constraining the phase curve if measured over a 
large enough range of phase angles, and hence over 
long enough observation arcs (C. Hergenrother, priv. 
comm.). The photometric S/N constraint is more strin¬ 
gent than the minimal one for astrometry, and thus as¬ 
trometric measurements could be made “for free” from 
these observations. Or, put another way, useful phase 
curve photometry could come out of the data taken for 
astrometry if the field were large enough to contain 
some of the growing number of faint photometric stan¬ 
dard stars. Given that there are already sufficient stars 
in a 30 arcmin field of view to permit photometry at 
the few percent level (Pickles and Depagne, 2010) and 
that the AAVSO Photometric All-Sky Surveyis still 
in progress, a 20 arcmin field like that suggested for as¬ 
trometry will prove sufficient. The photometric standard 
surveys have a faint limit of V ~ 16-17. While expo¬ 
sures targeting the magnitudes expected for the faintest 
NBA discoveries might have these faint standards satu¬ 
rated, interweaving quick exposures would permit pho¬ 
tometric calibration. 

From Table 1 we conclude that a 2 m telescope could 
be used for phase curve measurements for most H <22 
reapparitions as well as most new discoveries out to a 
couple of weeks. A 4 m could extend coverage to a 
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month for the mean new discovery at an hour per ex¬ 
posure (but keep in mind the complications from light 
curve modulation; see §8.5). This might not be justi¬ 
fiable based on the phase curve knowledge alone but it 
would make sense to plan the long-arc astrometric fol¬ 
lowup observations to have longer exposure times (ei¬ 
ther in single exposures or by coadding) than strictly 
necessary for the astrometry itself, to permit the phase 
curve measurements at large phase angles. 

9. Conclusions 

We have used properties of the over 10,000 known 
NBAs and their discovery circumstances to show that 
rapid follow-up after discovery is essential for closing 
the wide, and increasing, gap between discovery and 
characterization on a decade timescale, particularly as 
an increasing number of new discoveries will be of 
smaller {H > 22) NBAs that are projected to be lOOx 
fainter on subsequent apparitions. 

Spectral follow-up that keeps pace with discoveries 
could be accomplished with a dedicated 2 m telescope 
observing NBAs within days of discovery. Bong-arc 
(months) astrometry would require a 4-r m telescope, 
and the same observations could be used to measure 
phase curves and colors that would help to constrain 
albedos. Bight curves for the bigger (and likely more 
slowly-rotating) new NBAs could be collected shortly 
after discovery with a 2 m telescope, but the cadence re¬ 
quired to measure fast rotations likely for smaller new 
NBAs would necessitate rapid follow-up with a 4 m 
telescope. 

Already-known large {H < 22) NBAs making return 
apparitions are predicted to reach brightnesses similar 
to their discovery magnitudes (generally brighter than 
the discoveries projected for new surveys, which will 
thus likely measure astrometry for them by default). For 
this subset of objects, targeted follow-up spectroscopy 
and photometry can be done with 2 m telescopes. For 
the brightest of these, 1 m-class amateur telescopes can 
contribute light curves and astrometry. 

If means are insufficient to fund the necessary tele¬ 
scopes, which should be prioritized? Spectral classifi¬ 
cation yields the highest scientific return as it provides 
composition, which also yields an albedo estimate, and 
thus a more stringent constraint on an object’s size. 
Bight curves may provide structural information or re¬ 
veal a binary system. 

Bastly, astrometric measurements are of value if one 
desires to perform targeted observations at a subsequent 
apparition after discovery and thus requires a precise or¬ 
bit. It should be noted that orbits obtained within the 
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first few days or weeks of an NEA’s discovery are gen¬ 
erally good enough to securely establish the shape and 
location of the orbit in space (Lee et aL, 2015), while 
further observations, especially at future aparitions, will 
do most to improve the accuracy of locating the asteroid 
along its orbit, and thus on the sky. The level of accu¬ 
racy achieved for the initial orbit is sufficient for studies 
of the shapes and distribution of asteroid orbits and their 
evolution over time; it is also sufficient to determine if it 
intersects with, or is close enough to. Earth’s orbit that 
an asteroid with that orbit could be at risk of impact with 
our planet. 

Because the majority of new NBAs are not going to 
pose any danger of short- or mid-term impact, it may 
be argued that a more precise orbit at time of discovery 
is unnecessary because an NBA will be serendipitously 
reobserved by the surveys at some future apparition(s) 
and linked back by the Minor Planet Center to its dis¬ 
covery apparition. In this case, improving the accuracy 
of an orbit becomes a game of wait-n’-see, which ar¬ 
gues strongly for characterization to take place during 
the discovery apparition as there is no guarantee a new 
NBA will be found at a later apparition with a targeted 
observation. This strategy, however, is unacceptable for 
mission planning, where the location of a target aster¬ 
oid would need to be known with high accuracy years 
ahead of launch. Commercial ventures planning to send 
probes to many NBAs of varying sizes (Elvis and Esty, 
2014) would benefit greatly from more accurate orbits 
being determined during the discovery apparition if they 
intend to be able to act quickly on newly-discovered 
profitable asteroids. 

We have shown that a small group of dedicated tele¬ 
scopes (one 4 m and two 2 m, with help from smaller 
amateur contributions) could characterize or constrain 
the spectral class, albedo and spin period of most new 
NEA discoveries while also working through the back¬ 
log of known NBAs. There are compelling scientific 
reasons to do, which also have implications for plane¬ 
tary defence: Understanding the composition and struc¬ 
tural make up of the smaller {H < 22) NBAs will 
provide knowledge of the size group of asteroids most 
likely to impact Earth in the near future, while testing 
current theories of how the NEA population is fed by 
the Main Belt. 
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